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ABSTRACT 

High-resolution N-body simulations are used to examine the power spectrum dependence of the concentration 
of galaxy-sized dark matter halos. It is found that dark halo concentrations depend on the amplitude of mass 
fluctuations as well as on the ratio of power between small and virial mass scales. This finding is consistent with 
the original results of Navarro, Frenk & White (NFW), and allows their model to be extended to include power 
spectra substantially different from Cold Dark Matter (CDM). In particular, the single-parameter model presented 
here fits the concentration dependence on halo mass for truncated power spectra, such as those expected in the 
warm dark matter (WDM) scenario, and predicts a stronger redshift dependence for the concentration of CDM 
halos than proposed by NFW. The latter conclusion confirms recent suggestions by Bullock et al., although this 
new modeling differs from theirs in detail. These findings imply that observational limits on the concentration, 
such as those provided by estimates of the dark matter content within individual galaxies, may be used to constrain 
the amplitude of mass fluctuations on galactic and subgalactic scales. The constraints on ACDM models posed 
by the dark mass within the solar circle in the Milky Way and by the zero-point of the Tully-Fisher relation are 
revisited, with the result that neither dataset is clearly incompatible with the 'concordance' (f2o = 0.3, Ao = 0.7, 
us = 0.9) ACDM cosmogony. This conclusion differs from that reached recently by Navarro & Steinmetz, a 
disagreement that can be traced to inconsistencies in the normalization of the ACDM power spectrum used in that 
work. 

Subject headings: cosmology — dark matter — galaxies: formation — galaxies: structure — 



1. INTRODUCTION 

Due to their large density, the central regions of dark matter 
halos, where galaxies form according to the current paradigm 
of structure formation, hold important astrophysical clues to the 
nature of dark matter. This is why many studies have attempted 
to constrain dark matter models on the basis of clues to the dark 
mass distribution gained from detailed studies of the dynamics 
of gas and stars in individual galaxies. The most straightfor- 
ward method compares dark mass distributions derived from 
rotation curves of disk galaxies with detailed predictions of N- 
body simulations (Frenk et al. 1988; Flores et al. 1993; Flores 
& Primack 1994; Moore 1994; Moore et al. 199b), although 
similar insight can be gained by inspecting the high-order mo- 
ments of the stellar velocity distribution in spheroid-dominated 
systems (Carollo et al. 1995; Rix et al. 1997; Gerhard et al. 
1998; Cretton et al. 2000; Kronawitter et al. 2000). 

Despite the simplicity of the rotation-curve method and the 
numerous studies reported in the literature to date (see, e.g., 
Swaters 1999 for a comprehensive list of references), there is 
still no broad consensus regarding the detailed distribution of 
dark matter in disk galaxies, a situation that reflects the dif- 
ficulties associated with obtaining accurate circular velocities 
over a large dynamic range in radius, as well as with accounting 
for the contribution of the baryonic component to the rotation 
curve and for the uncertain response of the dark material to the 
assembly of the galaxy. For example, while constant-density 
'cores' in the dark mass distribution appeared at first to be nec- 
essary to explain the rotation curves of low surface brightness 
(LSB) dwarf galaxies (Flores & Primack 1994, Moore 1994, 
McGaugh & de Blok 1998), the persuasiveness of the obser- 
vational evidence for these cores has recently been called into 



question by careful reanalysis of the observational datasets (van 
den Bosch et al. 2000; Swaters, Madore & Trewhella 2000; van 
den Bosch & Swaters 2000). 

At the same time, there is also considerable uncertainty in 
theoretical predictions of the dark mass distribution at radii as 
small as those probed by the rotation curve data. Most workers 
agree that Cold Dark Matter (CDM) halos have density profiles 
that diverge near the middle (a result that would be at odds with 
the alleged cores of LSB dwarfs), but there is still controversy 
as to the exact asymptotic behavior of the density near r = 0. 
The work of Navarro, Frenk & White (1996, 1997) suggested 
that the central density may diverge as fast as r , but subse- 
quent work has argued both for steeper (e.g. r" 1 4 in Moore et 
al. 1998) and shallower profiles (e.g., r" 7 in Kravtsov et al. 
1998, although it should be noted that the authors have appar- 
ently now retracted this result, see Klypin et al. 2000). Each 
of these models predicts, of course, quite different dark matter 
contributions to disk galaxy rotation curves, making it difficult 
to provide a sound interpretation of the observational data. 

In other words, even if observations could constrain beyond 
the dark mass distribution near the middle of disk galaxies, then 
there would still be no consensus on the exact significance of 
that finding for dark matter models. The reasons for the dis- 
agreements in the theoretical predictions are still being investi- 
gated, but in all probability they reflect the inherent difficulties 
associated with simulating accurately and reliably the dynami- 
cal behavior within individual galaxies, where the density con- 
trast exceeds 10 6 . Particles inhabiting these regions go about 
their orbits thousands of times during a Hubble time, making 
numerical results highly vulnerable to insidious systematic ar- 
tifacts associated with the choice of integrator, time-stepping, 
and gravitational softening. Unfortunately a full account of the 
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dependence of the innermost density profiles of CDM halos on 
such numerical parameters is still lacking, but the indication is 
that it will require extreme care and a concerted numerical ef- 
fort on massively parallel computers to be able to characterize 
unequivocally the behavior of the dark matter density profile 
within the regions probed by rotation curve data. 

Given the intrinsic difficulty in providing robust theoretical 
predictions for the shape of the inner density profiles and the 
unsettled status of the interpretation of current rotation curve 
datasets, it is important to identify alternative observational 
and theoretical comparison criteria that are less sensitive to nu- 
merical and observational shortcomings. Navarro & Steinmetz 
(2000a,b, hereafter NS00a,b) have recently argued that one pos- 
sible choice is to use the total dark matter content within the 
main body of individual spiral galaxies. 

The typical radii involved are of order ~ 10 kpc for a bright 
spiral, which corresponds to about 3-5% of the virial radii. 
These regions are much less affected by numerical resolution 
issues than the ^kpc regions probed by rotation curves. Also, 
by focusing on the total dark mass within this radius rather 
than on its detailed radial distribution, both observational and 
theoretical estimates are presumably more reliable. For exam- 
ple, as discussed by NSOOa, there are strict upper limits on the 
dark mass enclosed within the solar circle in the Milky Way 
from detailed models of Galactic dynamics (Dehnen & Binney 
1998; Gerhard 2000). Such a constraint can be extended to 
other spiral galaxies by examining the zero-point of the Tully- 
Fisher (TF) relation. Indeed, provided that stellar mass-to-light 
ratios and exponential scalelengths can be estimated reliably, 
the TF relation allows for direct estimates of the dark mass 
within a couple of exponential scalelengths from the middle of 
the galaxy. 

NSOOa applied these constraints to a number of halos simu- 
lated within the ACDM scenario, and concluded that the dark 
mass in ACDM halos is too centrally concentrated to be consis- 
tent with observations. This result added to an uncomfortably 
long list of concerns regarding the viability of CDM on the scale 
of individual galaxies, including the survival of a large number 
of small mass halos within the virialized body of a parent halo 
(the 'substructure' problem, see Klypin et al. 1999; Moore et 
al. 1999a), as well as the evidence for constant density cores 
in dark halos alluded to above. Taken together, the evidence 
appeared to warrant a radical revision of one or more of the 
premises of the CDM paradigm, and there has been no shortage 
of proposals: self-interacting dark matter (Spergel & Steinhardt 
2000), warm dark matter (Hogan & Dalcanton 2000), fluid 
dark matter (Peebles 2000), fuzzy dark matter (Hu, Barkana 
& Gruzinov 2000), etc, all aim to provide a model that behaves 
like CDM on large scales but with reduced substructure and 
'concentration' on the scale of individual galactic halos. 

If the results of NS00a,b hold and ACDM halos are too con- 
centrated to be consistent with observations, then what changes 
are needed in order to reconcile the predictions of this sce- 
nario with observations? Are changes in the overall normal- 
ization of the power spectrum necessary, or does the shape of 
the ACDM spectrum require modification? Do small-scale cut- 
offs in the power spectrum (as expected in warm dark matter 
models) help? Or, in a more general sense, what is the rela- 
tionship between halo concentration and the power spectrum of 
initial density fluctuations? 

These are the questions addressed here through an extensive 
suite of N-body simulations. A description of the numerical 



simulations is given in Section 2, including details of the var- 
ious power spectra chosen for this study. Section 3 contains 
the main results regarding the concentration of dark matter ha- 
los and their dependence on the power spectrum, and Section 4 
uses these results to revisit the viability of the ACDM model re- 
garding the Milky Way and Tully-Fisher constraints. Section 5 
summarizes the main conclusions. 

2. NUMERICAL METHODS 
2.1. Cosmology and Power Spectra 

All of the simulations described here adopt the same cos- 
mological background model: a flat, cosmological constant- 
dominated universe with matter density parameter f^o = 0.3, 
Ao = 0.7, and Hubble parameter h = 0.65 1 . Two different power 
spectrum shapes have been considered. The first is the standard 
CDM spectrum, in the form given by Bardeen et al. (1986), 
which is fully characterized by erg, the present linear theory am- 
plitude of mass fluctuations in spheres of radius 8/i _1 Mpc, and 
by the value of the 'shape' parameter, T (Bardeen et al. 1986, 
Sugiyama 1995). 

The second power spectrum shape aims to mimic a warm 
dark matter (WDM) power spectrum: it is identical to the 
CDM spectrum on large scales but its power is reduced on 
scales smaller than that of a characteristic free-streaming mass, 
AvdmW = PcDM(k)exp(-kRf-(kRf) 2 )), where Rj is the co- 
moving free-streaming scale. Following Sommer-Larsen & 
Dolgov (2000) and Avila-Reese et al. (2000), a free-streaming 
wavenumber, kf, is defined as that where the WDM power spec- 
trum is half the value for CDM. This implies kf « 0.46 /Rf. The 

'The present value of Hubble's constant is parameterized by Hq = 
lOOftkms-'MpcT 1 . 
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FIG. 1. — The z = linear amplitude of mass fluctuations for the seven 
models investigated in this study. <r(M) is calculated using a top-hat real-space 
window function. Solid lines represent, from bottom to top, ACDM models 
So.9> Si. 2 and Si. 6. The short-dashed line corresponds to To. 5 and the dotted 
line to To.i . The two WDM models are shown with long-dashed lines; the top 
and bottom lines correspond to models W[ and Wg respectively. 



Eke, Navarro & Steinmetz 



3 



(1) 



free-streaming mass is defined as 

Mf = y PWDM I y 

with Xf = 2n/kf and pwdm being the density of WDM. Ex- 
pressing the free-streaming mass in terms of the free-streaming 
scale yields 

M f = 3.7 x lO 14 O o (^/A _1 Mpc) 3 r 1 M . (2) 

This approximation to the actual WDM cosmogony neglects the 
non-zero velocity dispersion of the warm dark matter particle 
candidates, but recent work indicates that this omission should 
have negligible consequences for the quantities of interest here 
(Avila-Reese et al. 2000; Bode, Ostriker & Turok 2000). On the 
other hand, one advantage of this approximation is that the only 
difference between the CDM and WDM runs is the small-scale 
behavior of the power spectrum, which implies that systematic 
trends of halo structure with power spectrum shape are easier 
to identify. 

Table 1 contains a list of the specific parameters chosen for 
the various models, and Figure 1 shows cr(M), the z = am- 
plitude of linear mass fluctuations in spheres of a given mass 
corresponding to each power spectrum. In total, seven differ- 
ent models were investigated; five ACDM models with dif- 
ferent parameter choices for erg and Y and two WDM mod- 
els with different free-streaming masses, Mf. Model S0.9 will 
be referred to hereafter as the 'fiducial' model, because it is 
roughly consistent with the local abundance of galaxy clusters 
(Eke et al. 1996) and with the amplitude of CMB fluctua- 
tions (Stompor, Gorski & Banday 1995; Liddle et al. 1996). 
While a value of T = 0.2 was adopted for this default model, it 
is worth noting that, according to the fit of Sugiyama (1995), 
T = fl h &xp(-n h -(h/0.5) 1/2 (fl h /no)), and therefore T w 0.16 
would be a more appropriate value for the high baryon density 
parameter, n b « 0.019/T 2 = 0.045 (for h = 0.65), advocated by 
Tytler et al. (2000). 

2.2. The Simulations 

For each model listed in Table 1, the AP 3 M code (Couch- 
man 1991) was used to evolve 128 3 dark matter particles in a 
32.5/i _1 Mpc cube from z = 24 to z = using 2000 equal steps 
in expansion factor. At z = 0, four halos with circular velocities 
between 180 and 230km s" 1 (similar to that of the Milky Way) 
were selected for resimulation from a list of halos identified 
by the spherical overdensity group-finding algorithm (Lacey & 
Cole 1994). Unless otherwise specified, halo circular veloci- 
ties, Va, are measured at the virial radius, r\\ the radius of 
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Power spectrum parameter choices. 
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a sphere containing a mean density A times the critical value. 
The parameter A depends on £1 and A according to (e.g. Eke, 
Navarro &Frenk 1998), 

'ft - 30 , ifA = 0; 
ft - 45 , ifO + A=l 



A(fl,A)= 178 



(3) 



and is w 100 at z = for the cosmology adopted here. In ad- 
dition to the circular velocity cuts, a criterion of relative isola- 
tion was also enforced, so that halos considered for resimula- 
tion were restricted to those without neighbors more massive 
than 2.7 x 10 11 /T 1 M Q within 1.5/! _1 Mpc. This selection cri- 
terion increases the likelihood that the selected halos are close 
to equilibrium, simplifying the interpretation of the results. Be- 
sides the four 'Milky Way' halos selected for each cosmogony, 
further halos extending to circular velocities of order 100 kms -1 
were also selected for resimulation in the fiducial S0.9 model 
and the WDM models. 

The resimulations were performed using a multiple time-step 
N-body code based on the algorithm described by Navarro & 
White (1993), modified to take advantage of the GRAPE3 hard- 
ware (Sugimoto et al. 1990). Particles were allowed to take 
up to 86,000 time-steps during their evolution from the start- 
ing redshift of 50 to z = 0. Each halo has between 35,000 and 
85,000 particles within the virial radius at the final time. A 
Plummer gravitational softening of e = 0.6 kpc was used in all 
resimulations of 'Milky Way' halos. The extra resimulated ha- 
los with Va < 160kms _1 were run using e = 0.4 kpc. A few 
simulations were rerun varying the numbers of particles, and 
indicate that this numerical setup is appropriate for making re- 
liable measurements of the total mass within 5-10 kpc. 

3. POWER SPECTRUM AND HALO CONCENTRATION 

3.1. CDM and WDM density profiles 

Figure 2 shows the density profiles at z = corresponding 
to the cosmologies listed in Table 1 . Each profile is an average 
over the four 'Milky Way' halos (i.e., halos with Va in the range 
180-230 km/s). The mean profile for each model is shown, 
together with fits of the form proposed by Navarro, Frenk & 
White (1996, 1997, hereafter NFW), 

^ = «E (4 ) 

Pent (r/r s )(l + r/r s ) 2 ' 

where p cr ; t = 3H 2 /8irG is the critical density for closure, 5 C is 
a characteristic density contrast, and r s is a scale radius that 
corresponds to the region where the logarithmic slope of the 
density equals the isothermal value, dln(p)/dln(r) = -2. 

The main point to note here is that the NFW fitting formula 
works quite well for ACDM halos in the radial range 0.1-10 r s , 
in agreement with the results of NFW. This fitting formula also 
reproduces the density profiles of WDM halos, even for mass 
scales well below the free-streaming mass, Mf. This result has 
been noted before (Huss, Jain & Steinmetz 1999; Avila-Reese 
et al. 2000; Bode et al. 2000), and allows the characterization 
of each halo by two simple parameters: the mass inside ta (the 
virial mass Ma) and the 'concentration' ca = fA/fs- The con- 
centration is directly related to the NFW characteristic density 
contrast by 

S -± c l ( 5 ) 

c 3 [ln(l+c A )-c A /(l+CA)]' 

so that either parameter describes fully the density structure of 
a halo of a given mass. In what follows, ca will be adopted 



4 



Power Spectrum and Dark Halo Concentration 



except in the comparison with the results of NFW, where 6 C 
will be used. This is motivated by the fact that NFW adopted 
A = 200 in their work, whereas the more general A definition of 
eq. 3 is adopted here. Note that S c is independent of A, but that 
concentration is not, so that one should be careful when com- 
paring concentration values quoted by different authors. For 
the model considered here, A « 100 at z = 0, and ca ~ 1 -3 c 2 oo- 
Note that, although the resolution of these simulations is good 
enough to measure concentrations in a robust manner, it is not 
adequate to address the ongoing controversy regarding the in- 
nermost slope of the density profile. 

3.2. The Mass Dependence of Halo Concentration 

Figure 3 shows the concentrations measured in the simula- 
tions at z = 0, as a function of the virial mass of each halo. 
Different symbols correspond to different cosmogonies, as de- 
scribed in the caption to Figure 2. The top panel corresponds 
to the three ACDM models, Si.6, Si. 2, and S0.9, from top to bot- 
tom, respectively. The middle panel shows models To. 5 and 
To.i, while the bottom panel presents results corresponding to 
the warm dark matter models W\ and W$. Results for the fidu- 
cial model Sq.9 are repeated in all panels. 

There are a few things to note in this figure. Firstly, ACDM 
concentrations increase with increasing a& and decrease with 
increasing mass. These trends are consistent with those re- 
ported by NFW on the basis of lower resolution simulations, 
and support NFW's interpretation that the concentration, or 
equivalently, the characteristic density of a halo, reflects the 
mean density of the universe at a suitably defined collapse time. 
Collapse redshifts increase for higher values of the normaliza- 
tion parameter erg, and are higher for low mass systems, re- 
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FIG . 2. — Density profiles of dark matter halos formed in the different cosmo- 
logical models. Different symbols correspond to different models, as follows: 
filled squares (S0.9), stars (S1.6), open squares (Si. 2), open circles CTo.i), filled 
circles (To.s), filled triangles (Wi), and open triangles (W%). Profiles shown are 
averages over the four 'Milky Way' halos (Va ~ 200 km/s) resimulated for 
each cosmogony. A second average halo is also shown in the lower right panel 
for models S0.9, W\ and Wg corresponding to Va ~ 100 km/s. The curves show 
NFW profiles fitted to the average profiles. Line types are as in Figure 1. 
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FIG. 3. — Concentration as a function of halo mass for all models, split into 
three panels for ease of presentation. Lines in each panel are fits to the data 
using the model described in §3.3, using the same value of C CT = 28 (the only 
free parameter in the modeling) for all models. The top panel shows, from 
top to bottom, ACDM models S\,e (starred symbols), S1.2 (open squares), and 
the fiducial model S0.9 (filled squares). The fiducial model is repeated in all 
panels for comparison. The middle panel shows models To. 5 (filled circles), 
and To.! (open circles). The bottom panel shows W\ (filled triangles) and Wg 
(open triangles). All line types are as in Figure 1. 
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fleeting the hierarchical development of structure in CDM uni- 
verses. 

Secondly, ACDM concentrations depend very weakly on 
mass for the range considered here; changing by only about 
50% over two decades in mass for model S0.9. Figure 4 shows 
the simulation results presented by NFW for a variety of differ- 
ent cosmological models and dark matter power spectra, P(k). 
The weak dependence on mass for the CDM models is surpris- 
ing when compared with the stronger trends observed for the 
power-law power spectra simulations, labeled with the spectral 
index n, where P(k) oc k n . As n becomes more negative, the 
concentration depends more weakly upon mass. This is to be 
expected, since, as pointed out by NFW, the scaling between <5 C 
and Ma found in their numerical simulation is S c oc 
the same that links the characteristic non-linear mass M*(z) and 
the mean cosmic density at redshift z? However, as can be 
readily seen in Figure 4, the S c -M dependence found for CDM 
models is actually much weaker than expected for n ~ -1.5, 
the 'effective' CDM spectral index on the mass scales probed 
by the NFW simulations. A more negative spectral index seems 
necessary to explain the CDM results. This led NFW to postu- 
late that it is the amplitude of fluctuations on mass scales much 
smaller than the virial mass that determine the concentration. 
Consequently they introduced a (rather arbitrary) parameter of 
order <; 1 % (see parameter / in equation 7 below) in their mod- 
eling, in order to shift the mass scale under consideration and 
reproduce the numerical results. This is a rather unsatisfactory 
aspect of their modeling that lacks clear interpretation. 

Finally, further clues can be gleaned from the concentration 
of To.5 and T0.1 halos. To.i concentrations are lower than S0.9, 
which is not surprising given that the amplitude of mass fluc- 
tuations is significantly lower on galactic scales (Figure 1). On 
the other hand, Tq, 5 concentrations are as high as S1.6, although 
<j(M) is in this case lower than Si. 6 on galaxy-mass scales (Fig- 
ure 1). This again hints that the amplitude on virial mass scales 
is a poor predictor of the concentration. These hints are con- 
firmed by the results of WDM model W%, which shows a clear 
reversal of the concentration versus mass trend on scales below 
a few times the free-streaming mass Mj. W$, concentrations de- 
crease with decreasing mass despite the fact that WDM a(M) 
increases towards low masses before saturating at M <C Mf 
(Figure 1). 

3.3. A model for the power spectrum dependence of the 
concentration 

Although the model proposed by NFW captures many of the 
qualitative trends shown in Figure 3 it suffers from two main 
shortcomings: (i) it introduces two arbitrary parameters whose 
interpretation remains unclear, and (ii) it predicts a redshift de- 
pendence for the concentration that is weaker than found in re- 
cent numerical simulations (Bullock et al. 2000). Bullock et al. 
propose an alternative prescription, also with two free parame- 
ters, that results in improved agreement between the predicted 
redshift dependence of concentrations and the results of the nu- 
merical simulations. Their model follows NFW in associating a 
halo's characteristic density with the average background den- 
sity at collapse time, but differs from NFW in the definition of 
characteristic density and collapse time. 

More specifically, NFW take the characteristic density of a 

2 The characteristic clustering mass M* is defined so that <r(M*)D(z) = <5 cr j t 
(= 1.686 for CI = 1, consult Lacey and Cole (1993) and Eke, Cole & Frenk 
(1996) for other values of fi and A). 



halo to be 5 C (see eq. 5), and use a constant of proportionality, 
C, to relate this to the mean background density at the collapse 
redshift, z c , according to 

8 c = Cn Q (l+z c ) 3 . (6) 

The collapse time is defined as that when, according to the Press 
& Schechter approach (Press & Schechter 1974; Lacey & Cole 
1993) , half the virial mass of the halo was first contained in 
progenitors more massive than a fraction / of the final mass. 
This implies that 



erfc < 



<5critfc)-<WZo) 

{ ^2[a 2 (fM)-a 2 (M)] 



(7) 



where z denotes the redshift at which the halo is identified, 
^crit(z) = ^crit(O) / D(z) is the spherical top-hat model linear over- 
density threshold and D(z) represents the linear theory growth 
factor. (D(z) = (1 +z)~ l if fio = 1, Ao = is conventionally nor- 
malized to unity at z = 0; formulae for other values of f2o and 
Ao can be found in Peebles 1980.) For their simulations, NFW 
found a good fit by adjusting the two free parameters to be 
C= 3000 and/ = 0.01. 

Bullock et al., on the other hand, choose the characteristic 
density, p s , to be such that 

4-7T 

M A =— r]p s , (8) 



and specify the collapse redshift, z c , solely in terms of a(M), so 
that 

D(z c )o-(FM A ) = 1.686, (9) 

where F = 0.01. Their second free parameter, K, relates the 
characteristic density to the background density via 

p s = K 3 A(Zo)p cr it(Zc) 

and feeds through into the concentration as 



(10) 



c A =K 



l+Zo 



(ID 



K = 4 provides a good fit to their ACDM simulation results. 

This model, like that of NFW, suffers from the introduction 
of two arbitrary parameters (F and K) whose interpretation re- 
mains unclear. Furthermore, the definition of collapse epoch 
given in equation 9 implies that for a truncated power spectrum 
such as WDM, halo concentrations will still increase monotoni- 
cally with decreasing mass, approaching a constant at M <C M/. 
This is at odds with the results presented in the previous section 
(see also Bode et al. 2000), which show that WDM halo con- 
centrations decrease on mass scales below a few times the free- 
streaming mass. These results strongly suggest that it is not 
only the amplitude of the power spectrum, but also its shape, 
that determine the concentration of dark matter halos. In par- 
ticular, only a modeling that includes such shape dependence 
will be able to reproduce the somewhat counterintuitive depen- 
dence of concentration on halo mass found for truncated power 
spectra such as W& (see Figure 3). 

After some experimentation, a simple model has been pro- 
duced that matches the mass dependence of halo concentra- 
tions for the simulations presented here. Furthermore, the same 
model also fits all of the original NFW results, whilst modifying 
the redshift dependence of concentrations so that they are com- 
patible with the recent results of Bullock et al. The new model 
has a single free parameter and is of more general applicabil- 
ity, since it can be applied to truncated power spectra, where 
Bullock et al.'s prescription fails. This model also removes the 
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FIG. 4. — Mass dependence of halo characteristic densities, as reported by Navarro, Frenk & White (1997) (solid circles), compared with the results of the model 
described in §3.3. Constant C a in eq. 13 has been chosen in each case so as to provide a good fit to the simulation results at M 2 qq ~ M* . Model SCDM corresponds 
to the former 'standard' CDM (H = 1, erg = 0.63, T = 0.5). Model CDMA has tt = 0.25, A = 0.75, T = 0. 19, and ct 8 = 1 .3. The rest of the panels correspond to 
power-law power spectra, P(k) oc V ; the value of n and f2o is listed in each panel. Masses are normalized to the characteristic clustering mass M„ , defined so that 
o-(M«) = <5 crl , (=1.686 for tt = 1). This corresponds to M„ = 1.6 X lO 13 /r'M0 for SCDM and M« = 4.1 X I0 l3 h- [ M Q for CDMA. Note that excellent fits can be 
obtained in all cases with similar values of the single free parameter C CT . 



Eke, Navarro & Steinmetz 



7 



need for the arbitrary, small mass fraction constant introduced 
by NFW and Bullock et al. (f=0.01 in equation 7 and F=0.01 
in equation 9) by postulating that the concentration of a halo is 
controlled by a combination of the amplitude and shape of the 
power spectrum. 

Consider the 'effective' amplitude of the power spectrum on 
scale M, defined by, 



a ell (A/) = a(M) ( -^M^-(M) 
dln(M) 



(12) 



This effective amplitude modulates a(M) so that, for WDM- 
like spectra, it decreases on mass scales smaller than a few 
times the free-streaming mass My. In broad terms, a given mass 
scale M collapses when D(z)<r(M) is at least unity. This time is 
controlled by the redshift evolution of the linear growth factor, 
D(z), appropriate for the cosmological model under considera- 
tion. Following this, the collapse redshift, z c , of a halo of mass 
M may be identified as 



D(zc)<T eS (M s ) = 



(13) 



where C a is a constant and M s is the mass contained within 
'"max = 2.17r s , the radius at which the circular velocity of an 
NFW halo reaches its maximum. The requirement for collapse 
that D(z c )cr(M s ) > 1 implies that C a > -dlncr/dlnM. For a 
power-law fluctuation spectrum with P(k) oc k n then this yields 
C a >6/(n+3). As in the models of NFW and Bullock etal., the 
mean density of the universe at the collapse redshift can then be 
used to calculate a characteristic density for the halo. Defining 
the characteristic density of the halo to be, as in Bullock et al. 
(see eq. 8), 



: A(z )p crit (z„)c A , 



(14) 



and setting this to equal the spherical collapse top-hat density 
at the collapse epoch, p sc , where 

Psc (z c ) = A(z c ) Pcrit (z e ) = A(z c /° ( * +Z , c)3 (15) 



^(Zc) 



yields, 



4 



A(z c ) 0(z o ) /1+z, 



A(z c ) 0(z c ) V l+z„ 



(16) 



Equations 13 and 16 describe the concentration of a halo of 
given mass, once the single free parameter in eq. 13, C a , has 
been specified. As the characteristic mass scale at which the 
effective amplitude of the power spectrum is evaluated depends 
on r s , and therefore on ca, equations 13 and 16 need to be 
solved iteratively to yield the combination of ca and z c 3 . 

This model reproduces, with roughly the same value of C CT , 
the results of the simulations presented here, all of the origi- 
nal results of NFW, as well as the redshift dependence advo- 
cated by Bullock et al. This is shown by the curves in Figure 
3, which show the result of applying the model, at z = 0, to 
the seven cosmogonies adopted in this study. Solid line types 
are used for the S models, short-dashed and dotted lines for 
To.s and To.i, respectively, while long-dashed lines are used for 
WDM. All of the curves use the same value for the propor- 
tionality constant in eq. 13, C a = 28. The model reproduces 
very well the trends with mass, normalization, and shape of the 
power spectrum seen here, including the counterintuitive trend 
towards lower concentrations seen in the low-mass halos. 

3 An algorithm to perform this calculation for CDM and WDM power spec- 
tra is available on request from the authors. 



Table 2 contains a list of concentrations for a I0 l2 h~ l M Q halo 
identified at z = in a variety of commonly studied cosmolog- 
ical models. This illustrates the interplay between Oo, Ao, eg 
and T. 

The model described above also reproduces the original re- 
sults of NFW quite well. This is shown in Figure 4, where the 
density contrast S c is plotted as a function of the mass enclosed 
within a A = 200 sphere, the parameters used by NFW. It is 
apparent from this figure that the model also reproduces the re- 
sults of the eight cosmogonies studied by NFW, including open 
models with O as low as 0. 1 , again with approximately a single 
value of the constant C a . 

3.4. The Redshift Dependence of Halo Concentration 

According to equations 13 and 16, the model predicts that, at 
fixed halo mass, in an Einstein-de Sitter cosmogony ca(M,z) oc 
(1+z). This relation agrees with the prediction of the model by 
Bullock et al. for the evolution of halo concentration. However, 
for low density universes the scaling with redshift is not the 
same as theirs, and it is therefore important to verify that it is 
still in good agreement with the numerical results. 

Figure 5 compares the predictions of all three different con- 
centration models with the results of the numerical simulations 
at z = and 2. The comparison includes all halos in the S0.9 
and Ws simulations with more than 2500 particles within r&. 
Concentrations labelled with an 'ENS' subscript in the top row 
correspond to the model presented here, 'B' to Bullock et al.'s 
(middle row), and 'NFW' for the NFW model predictions in the 
bottom panels. ENS concentrations use C a = 28 in eq. 13. 'B' 
concentrations use F = 0.01 and K = 4. NFW concentrations 
use C = 3000 and / = 0.01 . The typical halo mass range probed 
varies from lOH-lO'V M at z = to IO'^-IO'V'Mq at 
z = 2. 

The top two panels show that the model presented here pre- 
dicts a redshift dependence in good agreement with the simu- 
lation results, both for S .9 and Wg. The middle panels show 
that the Bullock et al. model also fits the results of the fidu- 
cial ACDM runs at z = and z = 2, but that their model fails to 
capture the mass dependence seen in the W% simulations. This 
illustrates the point that was made earlier that it is the effective 
normalization (equation 12), rather than simply a{M), that de- 
termines halo concentrations. The bottom row highlights the 
weak redshift dependence predicted by the NFW model com- 
pared with the simulation results, as noted by Bullock et al.; 
it slightly under-predicts the fiducial model concentrations at 
z = but over-predicts them at z = 2. 

In summary, the data in Figure 5 shows that the redshift evo- 
lution predicted by the model presented here is consistent with 
the simulation results. Even at z = 5, the highest redshift with 



Table 2 

Values of ca and C200 for a 10 12 /r' M halo identified at 

Z = IN FOUR COMMONLY STUDIED COSMOLOGICAL MODELS. 



n 


Ao 


h 


^8 


r 


CA 


c 200 


1 





0.5 


0.5 


0.5 


12.3 


11.7 


1 





0.5 


0.5 


0.2 


6.7 


6.4 


0.3 





0.65 


0.9 


0.2 


18.6 


15.2 


0.3 


0.7 


0.65 


0.9 


0.2 


12.0 


8.9 
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FIG. 5. — Comparison between concentrations measured at z = (left panels) and z = 2 (right panels), and the predictions of three different models. 'ENS' 
corresponds to the model presented in this paper (top panels), 'B' to that in Bullock et al (2000, middle panels), and 'NFW to concentrations computed using the 
procedure outlined in the Appendix of NFW (bottom panels). The filled-in squares correspond to So.9 halos, whereas open triangles correspond to Wg halos. The 
solid and dashed curves show the model predictions for So.9 and W$, respectively. 



Eke, Navarro & Steinmetz 



9 



simulation data in figure 1 1 of Bullock et al., the concentrations 
of 10 12 /i -1 M Q halos are c\ = 2.6 and 2.5 for the model in Sec- 
tion 3.3 and that of Bullock et al. respectively. Thus, it is not 
possible to discriminate between the slightly different redshift 
dependences of these two models. 

4. COMPARISON WITH OBSERVATIONS 
4. 1 . ACDM and the Dark Mass within the Solar Circle 



C8 V 2 20«dark < 0.8, (18) 

where V220 is the circular velocity, V200, of the Milky Way halo 
in units of 220 kms" 1 , and m dark is M dark /(4.3 x 10 10 M Q ). This 
suggests that, for T = 0.2, the circular velocity of the halo of the 
Milky Way should be somewhat less than 220km s" 1 , unless 
(78 < 0.9 or the Milky Way halo has an unusually low concen- 
tration for its mass. The fiducial S0.9 model may be reconciled 
with the Milky Way constraint if V200 <; 1 95 km s" 1 . A knock-on 
effect of moving the Milky Way into a smaller halo would be 
that any predictions of the galaxy luminosity function made by 
mapping mass into luminosity using the properties of the Milky 
Way would find that Milky Way-type galaxies were more abun- 
dant than before. As M oc V 3 , and assuming that the luminosity 
of a galaxy is proportional to its mass, it only takes a 20% de- 
crease in halo circular velocity to halve the luminosity. Thus, 
it remains to be seen whether assigning the luminosity of the 
Milky Way to the many halos with V200 ~ 195 kms" 1 is con- 
sistent with the luminosity function of bright spirals (see, e.g., 
Cole et al. 1994; Cole et al. 2000). 




1.8 2 2.2 2.4 



10 glo( V 20o/ kmS_1 ) 

FIG. 6. — Dark matter halo masses within the solar circle (8.5 kpc) com- 
pared with constraints derived from dynamical observations of the Milky Way 
(hatched region). A vertical arrow marks a circular velocity of 220 km s~' . The 
different symbols correspond to different cosmogonies, as specified in Figure 
2. Lines correspond to NFW profiles assuming concentrations given by the 
model in §3.3. 



Observations of star and gas kinematics provide well de- 
fined constraints on the dark matter content of the Milky Way 
within the solar circle, Rq = 8.5 kpc. As discussed by NSOOa 
(see that paper for full references), a simple upper limit on 
the dark mass within Rq, Md- d &(r < Rq), may be obtained by 
combining the observed circular velocity at the Sun's location, 
V c (Rq) = 220 kms -1 , with estimates for the total mass and ex- 
ponential scalelength of the Galactic disk (Mdi s k = 6 x 1O 1O M ^— - N 
and rdisk = 3.5 kpc, respectively). The result (note that there is © 
a typo in eq. 1 of NSOOa), ^ 

M d . dA (r < R ) £ 4.3 x 1O 1O M , (17) ^ 

constitutes an upper limit because the simple calculation de- ^ 



scribed above neglects two potentially important effects: (i) the 
contribution of the bulge, and (ii) any potential contraction that 

the dark halo may have experienced as a result of the assembly " * 

of the galaxy. CO 

The constraint expressed in eq. 17 is straightforward to com- W 

pare with the results of the numerical simulations described j_, 

here. This is done in Figure 6, where the dark mass within s — ^ 

8.5 kpc is plotted as a function of V200, the halo circular ve- p 

locity for A = 200 used in NSOOa, for all of the halos. The ^ 

top panel shows halos formed in the S .9> S1.2. an d S1.6 models, v — s 

the three different normalizations chosen for the ACDM see- 3 

nario. At V200 = 220km s" 1 (marked with an arrow in Figure 6), tjjfj 

A^dark('" < Ro) increases approximately in proportion to as. In O 

light of the modeling described above, this can be attributed to ' 1 
the higher average collapse times that result from the choice of 
higher normalizations. 

The lines in this panel correspond to the mass within 8.5 kpc 
predicted by the model described in §3.3. Combining this 
model with the constraint in eq. 17 it is possible to estimate the 
range of circular velocities allowed for the halo of the Milky 
Way as a function of the normalization parameter e^: 
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A lower bound on the mass of the Milky Way halo may be 
derived by requiring that the total baryonic mass of the Galaxy 
does not exceed the baryon mass within the virial radius of 
the halo. Assuming n b = 0.019 /T 2 w 0.045, the minimum 
halo mass corresponds to a circular velocity of ~ 105kms _1 . 
This is lower than the 130kms _I derived by NSOOa, because of 
the lower baryon fraction (f^ = 0.0125/T 2 ) and slightly higher 
Hubble constant (h = 0.7) adopted by those authors. 

From Figure 6, halos in the fiducial S0.9 ACDM model with 
circular velocities in the range (105, 190)kms _1 appear consis- 
tent with the Milky Way constraint. For erg = 1 .2 the range of ac- 
ceptable halo masses is narrower, and essentially no halo agrees 
with the observational constraints if erg = 1.6. This is reminis- 
cent, although less stringent, than the conclusion reached by 
NSOOa, who argued that erg = 1.14 ACDM halos were too con- 
centrated to be consistent with this constraint. However, a re- 
analysis of the NSOOa dataset reveals that because of inconsis- 
tencies in the normalization procedure for their ACDM simu- 
lations 4 , those authors had effectively normalized their power 
spectra to a% ~ 1.6 rather than erg = 1.14. After correcting for 
this error, the results in NS00a,b are consistent with those re- 
ported here. 

4.2. The zero-point of the Tully-Fisher relation 

As discussed by NS00a,b, the analysis of Section 4.1 can 
be extended to other spiral galaxies by examining the corre- 
lation between galaxy luminosity and the rotation speed of 
their gas and stars: the Tully-Fisher relation (Tully & Fisher 
1977). Provided that stellar mass-to-light ratios and exponential 
scalelengths can be estimated reliably, it is possible to evaluate 
the disk contribution to the circular velocity at 2.2 exponen- 
tial scalelengths (where the disk contribution peaks and optical 
Tully-Fisher velocities are typically measured) and derive con- 
straints on the total dark mass contained within this radius. 

The more concentrated a halo is, the faster a disk of given 
mass and radial scale must rotate to attain centrifugal equi- 
librium. Thus, as shown by NS00a,b, the zero-point of the 
Tully-Fisher relation provides a direct constraint on halo con- 
centrations. Although these authors conclude that erg = 1.14 
ACDM halos are too concentrated to be consistent with the 
I-band Tully-Fisher relation, as discussed in §4.1 their con- 
clusions were affected by an inconsistent normalization of the 
power spectrum. Given that the simulations of NS00a,b effec- 
tively probed a erg 1=3 1 .6 ACDM model and concentration de- 
pends strongly on erg, it is appropriate to revisit the issue and 
verify whether the fiducial S0.9 model is consistent with obser- 
vations. 

4.2.1. Gasdynamical simulations 

To this aim, a number of gasdynamical simulations including 
star formation and feedback have been run using GRAPESPH, 
a code that combines the hardware N-body integrator GRAPE 
with the smoothed-particle hydrodynamics (SPH) technique 
(Steinmetz 1996). The simulation setup and analysis are identi- 
cal to those described in Navarro & Steinmetz (1997) and Stein- 

4 This error originated in the fact that NSOOa used the transfer function pro- 
posed by Davis et al. (1985) to displace particles, while normalizing the power 
spectrum using the value at the Nyquist frequency of the original low-resolution 
simulation given by the CDM transfer function fit of Bardeen et al. (1986). At 
this small scale, the two fits give power spectrum values that differ by almost a 
factor of two, and this led to a systematic discrepancy between actual and in- 
tended normalizations. This error only affected the ACDM models of NS00a,b. 
All other models, including NFW's, are free from this problem. 



metz & Navarro (1999), and the reader is referred there for de- 
tails. In brief, the same initial conditions described above for 
the dark matter-only runs are used with the addition of gas, as- 
suming a value of = 0.0125/T 2 for the baryon density param- 
eter. Models with gas have typically 20,000 gas particles and 
the same number of dark matter particles. Up to 5, 000 of these 
end up in a galaxy at z = 0, resulting in lower resolution than the 
N-body simulations discussed in previous sections. Gas parti- 
cle masses range from 4.5 x 10 6 M Q to 2.5 x 10 8 M Q , depending 
on the model considered. 

Model galaxies are unmistakably identified in the runs as star 
and gas clumps with high density contrast. Only halos with 
more than 500 dark particles within the virial radius have been 
retained for analysis. The properties of the luminous compo- 
nent are computed within a radius, r ga i = 20(V2oo/220kms" 1 )/r 1 
This radius contains all of the baryonic material associated with 
the galaxy and is well outside the region compromised by nu- 
merical resolution effects. All of the rotation speeds are also 
computed at that radius. Note that r g . d \ exceeds the radii at which 
Tully-Fisher velocities are typically measured, but given the 
lower resolution of these simulations, rotation speeds at smaller 
radii are quite uncertain. This comparison therefore assumes 
that the circular velocity curves of actual disk galaxies remain 
approximately flat out to r ga \. 

4.2.2. The I-band TF relation 

Figure 7 compares the observed I-band Tully-Fisher relation 
(dots) with the numerical results for galaxies selected in the 
fiducial 50.9 ACDM model (filled squares) and in the Wi WDM 
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1 .6 
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M,-5 log(h) 

FIG. 7. — The I-band Tully-Fisher relation compared to the result of numer- 
ical simulations in the fiducial <r g = 0.9 ACDM (filled squares) and in the W[ 
warm dark matter (filled triangles) scenarios. Dots are a compilation of the 
data by Giovanelli et al (1997), Mathewson et al (1992), and Han & Mould 
(1992). The solid line is the best fit to the data advocated by Giovanelli et al. 
Horizontal 'error bars' in the simulation results span the range in luminosities 
derived from assuming a Scalo or a Salpeter IMF. Simulation circular veloc- 
ities are measured at r„. d \ = 20(V20o/220kms~')/r' kpc. Note that the slope, 
scatter and zero-point of the numerical TF relation are all in reasonable agree- 
ment with observation. 
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model. There is reasonable agreement between observation and 
simulations. The slope of the numerical TF relation is consis- 
tent with the observed value, and the scatter is much smaller 
(0.12 mag rms for So.9 and 0.10 mag rms for W\) than ob- 
served. These two conclusions are in agreement with the results 
ofNS00a,b. 

The main difference with that work is that now even the zero- 
point of the numerical relation appears to match reasonably 
well the observed value; the zero-point offset between simula- 
tions and observations is <~ 0.5 mag, compared with the 1 .5 mag 
offset reported by NS00a,b. The reason for the discrepancy can 
again be traced to the lower concentrations of Srj.9 halos com- 
pared with the results of NS00a,b. 5 Figure 7 also shows that 
there is little difference in the TF results obtained for W\ or S0.9, 
supporting the interpretation that the halo concentration is the 
main factor responsible for the zero-point of the numerical TF 
relation. 

The 0.5 mag difference between simulation and observation 
is not too worrying given that the simulated galaxies have col- 
ors that are slightly too red compared with their TF counter- 
parts. The average B-R color of the simulated galaxies is 1.2, 
with little dependence on luminosity. For comparison, the aver- 
age B-R in Courteau's (1997) sample is <~ 0.8. This suggests 
that star formation in the simulations occurs too early. Any 
modification to the feedback algorithm that remedies this will 
also tend to make the stellar population mix in the simulated 
galaxies brighter. If this correction can bring the stellar I-band 
mass-to-light ratios down from 2.5 to 1 .5, a value more in keep- 
ing with the results of Bell & de Jong (2000), then the 0.5 mag 
gap should be possible to bridge. 

In summary, it appears that if the I-band stellar mass-to-light 
ratio of TF galaxies is of order (M/L)j rts 1.5 then ACDM ha- 
los are consistent with the slope, scatter and zero-point of the 
I-band Tully-Fisher relation. Note however, that while halos 
formed in the fiducial ACDM scenario appear to have concen- 
trations consistent with observational constraints, other prob- 
lems associated with the assembly of disk galaxies through 
merging persist. In particular, the angular momentum (and 
size) of simulated disks is still quite below observed values, 
again suggesting that perhaps the feedback algorithm is not ef- 
fective enough at preventing the early collapse of baryons into 
protogalactic potential wells (Navarro & Steinmetz 1997). Ac- 
counting simultaneously for the luminosity, velocity, and an- 
gular momentum of spiral galaxies in these models remains a 
challenging problem for the ACDM cosmogony. 

5. SUMMARY AND CONCLUSIONS 

This paper contains the results from an extensive suite of nu- 
merical simulations which were aimed at understanding the re- 
lationship between the power spectrum of initial density fluc- 
tuations and the concentration of virialized dark matter halos. 
These simulations demonstrate that dark halo concentration de- 
pends both on the amplitude of mass fluctuations as well as on 
the shape of the power spectrum. A simple model that takes 
this into account by defining an effective amplitude as a(M) 
times the logarithmic derivative of a(M) with respect to mass 
on scales similar to the characteristic mass of the halo (i.e. that 
enclosed within the radius where the circular velocity peaks, 
'"max = 2.17r. s ) has been developed. This model reproduces the 
mass and redshift dependence of the concentration in all 7 cos- 

5 The normalization problem only affected the ACDM runs in those papers. 
All the results concerning the standard Q = 1 CDM model remain unchanged. 



mogonies investigated here, as well as in the 8 different cos- 
mogonies probed by NFW. It also extends the earlier models 
of NFW and Bullock et al. (2000) to power spectra very differ- 
ent from CDM, including truncated power spectra such as those 
appropriate for WDM. 

These findings are applied to the Milky Way, where obser- 
vational limits on the dark mater content within the solar cir- 
cle can be turned into constraints on the shape and normal- 
ization of the power spectrum. For the popular ACDM spec- 
trum, the Milky Way halo mass and the normalization of the 
power spectrum must satisfy the condition, erg V220 »*dark < 0.8, 
where V220 is the circular velocity of the halo (V200) m units of 
220 km s" 1 and m^k is the upper limit on the mass of dark mat- 
ter within 8.5 kpc of the middle of the Milky Way in units of 
4.3 x 1O 1O M . For erg = 0.9, the normalization favored from the 
abundance of galaxy clusters and CMB studies, this implies that 
the Milky Way halo has a circular velocity significantly smaller 
than the rotation speed at the solar circle, V200 < 195kms _1 . 
This finding may have significant impact on the luminosity 
function expected in this model, since 195kms _1 halos are 
much more abundant than their 220 km s -1 counterparts (see, 
e.g., Cole et al. 1994, 2000). Gasdynamical simulations includ- 
ing star formation and feedback also show that, because of their 
lower concentration relative to the NS00a,b study, ACDM ha- 
los are also roughly consistent with the zero-point of the I-band 
Tully-Fisher relation. The slope and scatter of this relation is 
also in good agreement with observed values. 

Halo concentrations in ACDM simulations are much lower 
than found by NS00a,b, who had argued that ACDM halos were 
too concentrated to be consistent with observations of the dy- 
namics of spiral galaxies. A reanalysis of their dataset reveals 
an inconsistency in the normalization of the power spectrum 
used in that work. Instead of the intended cr% = 1 . 14, their sim- 
ulations had an effective normalization of erg w 1.6. Once this 
correction is taken into account both studies yield consistent 
results. 

The set of simulations reported here thus identify and illus- 
trate the tight relation between power spectrum and halo con- 
centrations. The application of these results to the Milky Way 
and I-band Tully-Fisher relation lifts previous concerns and 
suggests that the concentration of erg = 0.9 ACDM halos is not 
clearly incompatible with observations. 
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